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ABSTRACT 

Context. Many solar filaments and prominences show short-lived horizontal threads lying parallel to the photosphere. 
Aims. In this work the possible link between Rayleigh-Taylor instabilities and thread lifetimes is investigated. 

Methods. This is done by calculating the eigenmodes of a thread modelled as a Cartesian slab under the presence of gravity. An 
analytical dispersion relation is derived using the incompressible assumption for the magnetohydrodynamic (MHD) perturbations. 
Results. The system allows a mode that is always stable, independently of the value of the Alfven speed in the thread. The character 
of this mode varies from being localised at the upper interface of the slab when the magnetic field is weak, to having a global nature 
and resembling the transverse kink mode when the magnetic field is strong. On the contrary, the slab model permits another mode 
that is unstable and localised at the lower interface when the magnetic field is weak. The growth rates of this mode can be very short, 
of the order of minutes for typical thread conditions. This Rayleigh-Taylor unstable mode becomes stable when the magnetic field is 
increased, and in the limit of strong magnetic field it is essentially a sausage magnetic mode. 

Conclusions. The gravity force might have a strong effect on the modes of oscillation of threads, depending on the value of the 
Alfven speed. In the case of threads in quiescent filaments, where the Alfven speed is presumably low, very short lifetimes are 
expected according to the slab model. In active region prominences, the stabilising effect of the magnetic tension might be enough to 
suppress the Rayleigh-Taylor instability for a wide range of wavelengths. 
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1. Introduction 

Threads are thought to be the building blocks of solar filaments 
and prominences. These fine-scale structures are readily seen in 
high quality images obtained using the Swedish Solar Telescope 
(SST), the Dutch Open Telescope (DOT), or the Solar Optical 
Telescope (SOT) onboard Hinode. Threads are a part of a mag- 
netic tube filled with cold plasma with densities much larger than 
the coronal surrounding, the density contrast being of the order 
100. Their observed widths are typically between 100 km and 
600 km, while th eir lengths are in the range 3500-28000 km (see 
' iLin et all 120051) . The total length of the magnetic tube where 
the thread is supposed to be embedded is much longer than the 
thread length, and in many cases is of the order of 10 5 km. Some 
threads are essentially horizontal, i.e., running parallel to the 
photosphe re, and are usually loc ated along the spine of the fila- 
ment (see Okamot o et all 120071 for a clear example). Inclined 
threads are mostly found in filament barbs, whereas vertical 
threads a re clearly observed in hedgerow pr ominences (see for 
example iBerger et all 120081: IChaeetall 120081) . 

Threads have short lifetimes, and in quiesce nt filament s they 
are in the range 9-20 min (see the review of lLinl |201 1). The 
mechanisms that produce the thread disappearance are usually 
thought to be related to mass flows and thermal instabil ities ( see 
for example, ICarbonell et ail l2004t ISoler et all 1201 ll 1201 2J in 
the context of prominences). However, other types of instabili- 
ties may play a role. Recently, Rv utova et all (1201 Oh have shown 
that a number of processes taking place in prominences during 
their evolution can be linked to fundamental fluid instabilities 
such as Rayleigh-Taylor and Kelvin-Helmholtz instabilities (see 



also lBereer et all [20101) . Using a single interface model between 
the prominence and corona, and the assumption of incompress- 
ibility for the perturbations, these authors have been able to at- 
tribute the appearance of plumes, spi kes and ripples to these in- 
stabilities. |HnnereFal] d2011l 2012|) have numerically investi- 
gated the nonlinear stability of the Kippenhahn-Shliiter promi- 
nence model to the Rayleigh-Taylor instability. 

they show 
Thompson & Schmieder ( 1991 ); 



An additional feature of threads is that 
oscillations. For example 



lYi & Engvoldl d 19911) : |YLeialJ dl99j found thatlome threads 
seem to oscillate independently whi le others appear to os- 
cillate togethe r (see t he reviews of Oliver & Ballesterl |2002; 
iBaneriee et all 120071: lArregui et ail |2012|) . More recent high 
spatial resolution and high cadence observations have revealed 
the presence of both propagating and standing waves in i ndivid - 
ual threads (s ee ILin et all I2007L l2009t lOkamoto et ail 120071) . 
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For example, ILin et all (120071) found oscillations with wave- 
lengths of the order of 3000 km and periods around 5 .4 min in 
quiescent filaments, while ILin et all (|2009) hav e reported peri- 
ods of 3.6 minutes in swaying filament threads. Okamoto et al. 
(2007) have detected waves in an active region prominence with 
wavelengths of at least 2 x 10 5 km with periods of 4 min. Note 
that the periods of oscillation are typically of the same order of 
magnitude as the lifetimes of threads in quiescent filaments. 

The theoretical interpretation of the reported oscillations 
in threads is usually done by means of magnetohydrody- 
namic (MHD) waves. A certain prominence model is as- 
sumed and the eigenmodes of oscillat i on are determined (see 
lOliver & Ballesterl .120021 lArregui et all I2012J for a detailed de- 
scription of the results using different configurations). In partic- 
ular, two basic mod els have been considered in the literature, th e 
Cartesian slab (see iJoarder eFall 1 1 997t IPiaz et all l200ll 120051) . 
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and the cylindrical tube (see for examp l e. | Di az et 



Dvmova & Rudermanl l200l ISoleri l20lol lArreguiet 
2011). MHD waves in cylindrical plasmas are characterised by 
two wavenumbers, the axial wavenumber and the azimuthal 
wavenumber. For typical threads the azimuthal wavenumber of 
the detected MHD waves is much larger than the longitudinal 
wavenumber. The advantage of the cylindrical tube is that it is 
more realistic than the Cartesian slab which is, in general, two 
dimensional and thus infinite in the third direction. However, the 
inclusion in the slab configuration of propagation in this third 
direction, equivalent to the azimuthal component in a cylindrical 
tube, mimics the effect of a three-dimensional structure regard- 
ing the MHD waves. 

The motivation of the present work is to study oscillations of 
horizontal threads under the presence of the gravitati onal force. 
This problem has been partially addressed by McEwan & Diaz 
d2007l) using a different approach. Here the thread is modelled 
using the slab configuration. The incompressible assumption is 
adopted for the perturbations which are allowed to be three- 
dimensional. The eigenmodes of the system are carefully anal- 
ysed, paying special attention to the regime where the system 
becomes unstable due to the Rayleigh-Taylor instability, driven 
by the gravity force. Our focus is on the possible relation be- 
tween this instability and the short lifetimes of threads. 



2. Single interface results 

It is well known that a single interface between fluids with dif- 
ferent densities, like the interface between threads and the sur- 
rounding corona, can be unstable in the presence of gravity when 
the upper fluid is heavier that the lower fluid, this is the Rayleigh- 
Taylor instability (Ravleigh, 1883; Tavlor, 1950). Let us assume 
that we have an interface in the z-direction. A constant magnetic 
field, Bo, is imposed to be along the jc-direction, and gravity, g, 
is acting across the magnetic field. The sketch of the configu- 
ration is represented in Figure [1] for a double interface, but at 
this point we just consider a single interface. We assume that the 
equilibrium density, po, is constant at each side of the interface 
(Po = Pp for z > 0, while po = p c for z < 0). Since the density is a 
piecewise constant magnitude the equilibrium gas pressure, po, 
has to change linearly with z due to the presence of the gravity 
force, and must be continuous at the interface to have hydrostatic 
balance. 



t/ 





+■ * 






g 


P c 








B 






+ 














J 

n 
\ 


1 
r 


/ 
/ 








P f 














L 


P c 





Fig. 1. Sketch of the slab model representing an idealised thread. 



We use the incompressible linearised ideal MHD equations 
for the perturbations. We seek for solutions with a dependence 
on t, x and y of the form e ,wt+( M+'* ) .y j where co is the frequency 
and k x and k y are the wavenumbers parallel and perpendicular 
to the magnetic field, respectively. The total wavenumber is k — 



Jk x + k 2 . The incompressibility constraint is 

dv, 
V ■ v = ik x v K + ikyVy + — — - 0. 
dz 

Under such assumption the continuity equation for the perturbed 
density reduces, in our configuration, to 



(1) 



icop — — V, 



dpo 

dz ' 



From the momentum equation we have 

itopov.x = -ik x p, 

iatpoVy = — ikyp H \ik x by - ik v b x ) , 

dp B I db x \ 

icop v z = — — + — ik x b z - —— - pg, 
dz p \ dz I 

while from the induction equation we obtain 

icob x — Boik x v x , 
icob y — Boik x Vy, 
icoby — Boik x v z . 



(2) 

(3) 

(4) 

(5) 



(6) 

(7) 
(8) 



Note that the equilibrium gas pressure, po, does not appear in the 
previous equations because of the incompressible assumption, 
but the perturbed pressure, p, is in the equations although we do 
not have an explicit equation for this magnitude. 

We want to obtain a single equation for v z , and the first step 
is to insert the corresponding components of the perturbed mag- 
netic field given by Eqs. ©-(O in Eqs. (@]l-((5]l and make use of 
Eq. (0. The following equations are obtained 

icopov y - -ikyp + (ik x Vy - ik y v x ) , (9) 

dp ^ B \k x [ dv x \ 1 dp 

iu>p v z = -— + ik x v z - — + —g— r v z . (10) 

dz p ti> \ dz J mo dz 

Eliminating the pressure perturbation from Eqs. (]3} and (0 we 
find that 



ik x v y — ik y v x — 0, 



(11) 



thus, the second term in the right hand side of Eq. (|9) is zero. 
Using Eqs. (0 and (0 together with Eq. (0 we find the follow- 
ing expression for the perturbed gas pressure 



-ipo 



co dv z 
¥~dl 



(12) 



In addition, combining Eqs. and ( TTTb we obtain 
k x dv z 

v * = l vTz- (13) 

Now inserting Eqs. (Tf2l and ( TOT l in Eq. (TTOb we find the final 
equation for v z 



p to 2 



dz \ dz 



dz 2 l V '- k [ P0+ c^^ U - {) 



(14) 



It is easy to see that for a piecewise density profile the eigen- 
functions have the following simple form 



v z (z) ■■ 



Ae kz , z<0, 
Be kz , z>0. 



(15) 
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The constants A and B are calculated requiring continuity of v z 
across the interface (meaning that A - B). There is an additional 
condition that comes from the integration of Eq. (fl4l i over an 
infinitesimal element of z that includes the interface (located at 
Z = 0). This condition reduces to 



-co pp+kr x — 

P 



dv z 
dz 



co 2 p t 



R2\ 



dv, 
dz 



-gk 2 (pp ~ Pc) v z (z = 0), 



(16) 



where we have used that po = p p for z > ("+" region) and 
Po = Pc for z < ("-" region). Inserting the derivatives of v, 
at each side of the interface in Eq. (fTol l leads to the well known 
dispe rsion relation for the i ncompressible interface problem (see 
also[Chandrasekhar] Il96ll) 



-gk 



Pp-Pc 



+ k 



2Bl 



Pp+Pc X ^{p v + Pc ) 



(17) 



From this equation it is clear that in the absence of magnetic 
field the system is always unstable (co 2 < 0) if p p > p c , i.e., 
when the heavier fluid is on top of the lighter fluid. The opposite 
configuration is always stable, i.e., co 2 > for p p < p c . The 
unstable mode can become stable if the second term of the right- 
hand side of Eq. ( fTTT i is larger than the absolute value of the first 
term. This term is simply the square of the kink frequency. It is 
interesting to note that Eq. (fTTT i simplifies to 



co - -gk + k x 



2^ 



(18) 



if the density at one side of the interface is much larger that the 
density at the other side (p p » p c , in the previous equation). This 
is a common characteristic in prominence threads. 

3. Slab model results 

The interface model is useful to understand the basics of the 
Rayleigh-Taylor instability. However, the main purpose of this 
work is to study threads, and these structures are better repre- 
sented using slab or cylindrical models. Here we assume a slab 
configuration. Across the field, i.e., in the z-direction, we as- 
sume again that density is constant. The equilibrium gas pressure 
is chosen to have hydrostatic equilibrium, and as in the interface 
model it changes linearly with z. The density takes the value p p 
inside the slab, representing a thread located between -zo and 
Zo, and p c outside the slab (see Figure [TJ. In this configuration 
the Alfven speed has a constant value inside and outside the 
thread. Notice that the thread is infinite in the y-direction. The 
model adopted here is similar to that studied in McEwan & Diaz 
(2007), where the main distinctive difference is the incompress- 
ible assumption for the MHD waves and specially the treatment 
of the boun dary conditions at the interf aces, given in our case 
by Eq. ( fT6b . Goedbloed & Poedts (2004) analysed a gravitating 
slab supported from below by a vacuum magnetic field but con- 
ducting w all boundary conditions were cons idered (see also the 
works of lHermans&Goossensiri987Ul989l) . 

The slab configuration is essentially the combination of two 
interfaces. Extending the single interface eigenfunction, given 
by Eq. tlSi . to the double interface problem, leads to 



(Ae k( - z+Zo \ z<-zo, 

v,(z) = | Be- k(z+Za) + Ce k(z - Z «\ z < |zol, 
[ De -k(z-zo) < Z>ZQ . 



(19) 



Continuity of v z at z ■ 

B = A-Ce- 2kz \ 



D 



C + Be 



-2kz 



+Zo imposes a relation between constants 

(20) 
(21) 



Since one of these constants can be chosen arbitrarily we assume 
that A - 1 . The constant D is written, using the equation for B, 
as a function of the constant C only. The next step is to apply the 
second boundary condition, given by Eq. (TToT l, at each interface 
(z = +Zo) using the appropriate values for the densities. These 
boundary conditions provide two equations for co 2 in terms of 
the constant C which is still unknown. The explicit form of this 
is calculated by combining these two equations. 

There are two solutions for C representing two different 
types of modes, hereafter denoted by the subscripts "+" and "-". 
The final dispersion relation is the following algebraic equation 



co 2 = 



-gk( Pp - Pc ) + kl 2 f(l-C ± e- 2kz «) 



Pp (l-2C ± e- 2kz °)+ Pc 



C ± = 



being 



where 
E + F 



*2*<"o 



gk 



1 + — |sinh2feo 
Pp) 



■k 2 v 2 



(22) 



(23) 



(24) 



.,2k 



i(*x P ) 2 



+8 2 k 2 



Pp) 



2 Pc 

sinh 2kzo H sinh 4kzo 

Pp 



1/2 



and 



G = 2smh2kz (k 2 y Ap +gk). 



(25) 



(26) 



In the previous expressions we have introduced the Alfven speed 
in the thread, va p = Bu/Jjlp p ~. It is clear that the slab dispersion 
relation given by Eq. (1221 has some similarities with the inter- 
face dispersion relation (Eq. (fTTt), but it is more involved due to 
the dependence of C on the different parameters. For a better un- 
derstanding of the solutions some limiting cases are investigated 
in the following subsections. 

Before we analyse the different cases we need to discuss the 
values of some of the parameters that appear in the dispersion 
relation. Although we are considering an infinite slab model in 
the y-direction we incorporate the three-dimensionality of the 
motion by choosing the appropriate k v . This is an important pa- 
rameter in our model. If we suppose that the thread is a mag- 
netic tube of circular cross-section then k v = m/R where m is 
the azimuthal wavenumber and R is the radius of the thread. For 
transverse motions m = 1 (since the system is periodic in the 
azimuthal direction m is restricted to be an integer). The equiv- 
alent mode in the slab has k y - 1 /zo (zo in the slab is equivalent 
to R in the cylinder) and corresponds to a transverse motion of 
the tube if k x <k k y . The longitudinal wavenumber is k x - 2n/A, 
being A the characteristic wavelength. It is clear that k x /k y is of 
order R/A. The typical radius of threads is R ~ 10 2 km, while 
the reported wavelengths are in the range A - 10 3 — 10 5 km. This 
means that the regime of interest is when R/A <K 1, i.e. in the 
thin tube approximation. 
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3.1. Purely magnetic case 

We consider first the case when gravity is absent and we anal- 
yse the properties of the eigenmodes. The system has symmet- 
ric or antisymmetric eigenfunctions, meaning that either B — C 
(for the "+" solution) or B — -C (for the "— " solution), and 
the eigenfunctions inside the slab are the hyperbolic functions 
"cosh" or "sinh", respectively. The symmetric mode, or kink, 
satisfies that D — A, while for the antisymmetric or sausage 
mode D - -A. In Figure [2] the eigenfunctions are plotted for 
a particular choice of parameters representative of thread con- 
ditions, R = 10 2 km (k y = 1/1?), A = 2 x 10 5 km (k x = 2n/A), 
p p - 100p c , and va p = l OOkms"" 1 . Contrary to the eigenmodes 
in a cy lindrical tub e (see ISpruitl[r9 82: Ed win & Robertsl 119831: 
Cally, 1986, 2003: iGoossens et all 120091) . the sausage mode in 
the slab model is trapped, instead of leaky. Additionally, the 
eigenfunction of the kink mode has a minimum at the centre 
of the slab instead of a maximum, but since the mode essen- 
tially produces a displacement of the whole slab , we sill refer to 
this mode as the transverse mode (see also Arregui e t al.l 120071 
for the spatial distribution of the eigenfunctions in the compres- 
sional problem). 



3 o.o 




Fig. 2. Sausage (top panel) and kink (bottom panel) eigenfunctions in 
the purely magnetic case. In this plot R = 10 2 km, A = 2 x 10 5 km, 
p p = 100p c , va p = 100 km s -1 andg = 0. The dashed lines represent the 
slab boundaries. 



It can be shown from Eqs. d23l-(l26li that in the case of zero 
gravity the constants are C± - 1/ fe~ 2fco ± 1 j, and the modes are 
always stable. The dispersion relation without gravity reduces to 



2 2 B\ 1 + coth(feo) 
to - k% 



for the kink mode and 
2 l2 B\ l+tanh(fe ) 



x fi Pp+Pc tanh(feo) ' 



(28) 



for the sausage mode. These equations agree with the results of 
Edwin & Roberts ( 1982) when k y — and in the incompressible 
limit (see their Eq. (12)). When kzo is large the hyperbolic tan- 
gent and cotangent tend to one, and we recover the dispersion 
relation for a single interface (see Eq. (ITTt ) when gravity is ab- 
sent, i.e., the frequency tends to the kink frequency. For a finite 
slab and when k x <k k y the frequency of these modes tends again 
to t he kink frequency (th is is true even for the compressible case, 
see lArregui et all [20071) . and this is the regime we are interested 
in. 



3.2. Purely gravitational case 

We now analyse the situation when gravity is present and Bq = 
0. For this case it can be shown that Eq. (l22l simplifies to the 
following dispersion relation 



co - ±- 



gk[p v -p c )sinh2fe 



-y(pp + Pc) sm h 2 2 feo + PpPc sinh4fe 



(29) 



x p p p + p c coth(feo) ' 



(27) 



This configuration has a stable mode always associated to C+, 
and an unstable solution linked to the C constant. The single 
non-magnetic interface dispersion relation, given by Eq. ([TBI 
without the magnetic term, is easily recovered from Eq. (|29l 
when p p » p c (the second term in the denominator can be ne- 
glected in front of the first term). If kzo » 1 the dispersion rela- 
tion tends also to the interface solution since the wavelength of 
the perturbation is much smaller than the half width of the slab. 
The spatial distribution of the corresponding eigenfunctions is 
rather different respect to the kink and sausage magnetic eigen- 
functions. It is possible to see that the stable mode is mainly 
localised at the upper interface since for this mode C - D (C+ 
tends to e~ 2kz ° according to Eqs. (I2"3l-(l2"6ii). while the unstable 
mode satisfies that B a A (C_ tends to 0), and it is therefore lo- 
cated at the lower interface (in this case the heavier fluid is on 
top of the lighter fluid and for this reason the mode is unstable). 
An example of the eigenfunctions is found in Figure [3] The de- 
gree of localisation of the eigenfunction around the interfaces is 
given by the penetration length, which according to Eq. (fT9l , can 
be defined as / = 1 /k. For the situation k x <k k v and k y - 1 /zo 
we have that / m zo (k ~ ky). Hence, for these parameters, slab 
eigenmodes due to gravity do not represent a coherent motion 
of the whole slab like in the kink mode, and are mostly surface 
waves associated to the individual interfaces. 

3.3. Full case 

When gravitational and magnetic forces are present at the same 
time we expect to find modes with mixed character. As in the in- 
terface problem, the interplay of these two restoring forces might 
cause the unstable mode to become stable if the magnetic field 
is strong enough, or if the wavelength is larger than a critical 
value. This behaviour is associated to the C solution in our no- 
tation. On the contrary, the stable mode found in the situation 
with gravity only, remains always stable when the magnetic field 
is different from zero. In this case the solution corresponds to the 
C+ constant. An example of this behaviour is found in Figure |4] 
where the square of the frequency calculated using the full dis- 
persion relation given by Eq. (l22l . is plotted as a function of 
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Fig. 3. Eigenfunction of the modes associated to the lower unstable 
mode (top panel) and upper stable mode (bottom panel) in the purely 
gravitational case. In this plot R = 10 2 km, A = 2 x 10 5 km, p p = 100p c 
and VAp = (g = 0.274 km s -2 ). The dashed lines represent the slab 
boundaries. 



the wavenumber. We identify the two modes, the stable solution 
(dashed line) and the initially unstable mode (continuous line) 
that becomes stable when the wavenumber is larger than a crit- 
ical value. We have used the same choice of parameters, mean- 
ing that we are in the regime k x <s fey. The solutions associated 
to the individual interfaces, calculated using Eq. (TPTl i. are also 
plotted in Figure |4] with thin lines. The agreement between the 
slab and interface solutions is quite good for small lc x , meaning 
that for the particular parameters we are using, the modes are, 
to all practical purposes, surface waves associated to the indi- 
vidual interfaces in this regime. Thus, for k x small there are no 
global kink-like transverse modes. Further confirmation of this 
behaviour is found in the eigenfunctions, represented in Figure 
[3] for va p = 100 km s _1 and A - 2x 10 4 km (see continuous 
lines). The surface waves are barely affected by the presence of 
the other interface, and the overall structure is quite similar to 
the one found in the non-magnetic case (compare with Figure 
0. 

The weakly coupled behaviour of the two interfaces mani- 
fested in Figure [5]is better understood by analysing the value of 
C in the dispersion relation. If we assume that p c <k p p the terms 
proportional to p c /p P can be neglected in Eqs. d24b and d25b . It is 
easy to see that with this simplification, and when the following 
condition is satisfied 



gk sinh 2fe » k\v 2 



Ap> 



the constants E and F tend to same value. This means that for 
the unstable mode the constant C goes to zero, and thus for 



Fig. 4. Square of the frequency as a function of k x . In this plot R = 
10 2 km, p p = 100p c , and v Ap = lOOkms -1 . Thick lines correspond 
to the full dispersion relation for the slab (Eq. 11221 ). while thin lines 
represent the solutions of the single interface (Eq. 1 1171 ), The dashed 
lines correspond to the C+ solution while the solid lines represent the 
C_ solution. The dotted vertical line represents the critical wavenumber 
calculated using Eq. ( 132b . 



this mode the slab dispersion relation tends to the interface dis- 
persion relation given by Eq. (fTST i. The equivalent result is found 
for the stable mode. Therefore, if the condition given by Eq. (l30l > 
is satisfied, the system is well described by the interface results. 
Note that Eq. (l30l does not mean that the role of the magnetic 
field is neglected in front of the gravitational term, the inter- 
face dispersion relation still contains the magnetic term. In this 
regime we can make further analytical progress assuming that 
k x <s k y , and reducing Eq. ([T"8T >. for the unstable mode, to 



R 



+2ky Ap , 



(31) 



where we have used that k y — 1 /zo — 1 //?. 

If we assume that the thread radius and the Alfven speed 
are known from observations, we can determine the critical 
wavenumber required to have a stable configuration (w 2 > 0) 



(32) 



In Figure |4] this critical wavenumber is represented by a verti- 
cal dotted line. Alternatively, if the wavenumber is known the 
Alfven speed required to have a stable configuration is 




K Ap 




(33) 



This equation might be useful for some applications, as we will 
show in section [4] 

The condition given by Eq. (f30b is not fully satisfied if, for 
example, we chose a wavelength of the order of A ~ 6 x 10 3 km, 
and we consider that the Alfven speed is va p = 100km s -1 , i.e. 
when we move to the right in the dispersion diagram shown in 
Figure[4] For these parameters the modes are no longer pure sur- 
face waves associated to the interface and have a mixed nature, 
as can be appreciated in Figure[5](see dotted lines). In fact in this 
figure the particular choice of parameters satisfies that 



(30) gfesinh2fco 



*Mp. 



(34) 



This condition means that the coupling between gravitational 
and magnetic forces is strong, and the resulting modes a mixture 
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Fig. 5. Eigenfunctions in the presence of gravity and magnetic field 
associated to the C_ solution (top panel) and the C+ solution (bottom 
panel). The continuous line correspond s to A = 2 x 10 4 km, the dotted 
line to A = 6 x 10 3 km, while for the dash-dotted curve A = 2 x 10 3 km. 
In this plot R = 10 2 km, p p = 100p c and v Ap = 100 km s _I . The dashed 
lines represent the slab boundaries. 



between interface and slab waves. The situation is intermediate 
between the purely gravitational and the purely magnetic modes. 
Finally, when 

gksmh2kz «ik 2 x v 2 Ar (35) 

the modes have a global nature and are very similar to the 
sausage and kink modes (see Figure [5] dash-dotted lines) since 
they are dominated by magnetic forces. This case is almost 
equivalent to the purely magnetic situation discussed in subsec- 
tion l3.ll and the modes are always stable. 

Figure [5] shows very clearly the transformation from sin- 
gle surface waves to global sausage and kink modes when the 
Alfven speed is fixed and the wavelength is allowed to change, 
but a completely equivalent behaviour is found if the wavelength 
is fixed and the Alfven speed varies. Note also that according to 
the dispersion diagram given by Figure 2] we find that the so- 
lution associated to the C constant changes from unstable to 
stable, and the eigenfunction evolves from an interface solution 
associated to the lower interface in the gravitationally dominated 
case to a sausage mode in the magnetically dominated situation. 
On the other hand, the solution associated to the upper interface, 
the C+ mode, is always stable and it eventually tends the global 
kink mode. 



3.4. Thread lifetimes 

From the previous results we conclude that in case of instability 
it will be triggered at the bottom of the thread and it is due to the 



C solution. The development of this instability will eventually 
produce significant changes in the whole thread structure. Here 
we want to calculate the time required to observe the effects on 
the upper interface. We know that the velocity of the lower inter- 
face is (in the linear regime) 



v e 



t/r 



(36) 



where t - 1 /\o>\ is the growth rate of the instability. If the inter- 
face moves upward then the vertical displacement is 



£ = v Te' /T + ft 



(37) 



We assume that at t — the interface is located at z = -Zo ( ac ~ 
cording to our slab model), this means that 



ft = - (zo + vqt) . 



(38) 



Now using Eq. ( l37l > we can calculate the time required to have a 
displacement of the interface from -zo to zo 



v r 



(39) 



This parameter provides a rough measure of the time to blend the 
lower and upper interfaces of the thread and, as a result of this 
process there is a mixture of the thread density with the much 
lower coronal density. The time scale t m can be interpreted as the 
disappearance time or lifetime of the thread due to the Rayleigh- 
Taylor instability. However, we have to bear in mind that the pre- 
vious calculation is based on the assumption of a linear regime, 
but in a real situation the nonlinear effects might increase the 
lifetimes due to the saturation of the instability. 

4. Application to observations of oscillating threads 

iLin et alJ (120071) found travelling waves in a quiescent fila- 
ment from the analysis of velocity Dopplergrams. The observed 
threads have a radius of the order of 100 km, and their typical 
lifetime is 9min. For the particular thread marked with the la- 
bel L2 in ILin et all d2007t) . the period is 5.4 min, the estimated 
wavelength of the propagating waves is 3000 km, and the phase 
speed is around 9.5 km s -1 . The amplitude of the oscillations is 
l-2kms _1 . 
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Fig. 6. Square of the frequency as a function of the Alfven speed in the 
thread. In this plot R= 10 2 km, A = 3 x 10 3 km and v^ p as 17.7 km s" 1 . 
The same notation as in the previous figures has been used. The long 
dashed line corresp onds to th e frequency (P = 5.4 min) derived from 
the observations of iLinet all ((2007). The dashed lines correspond to 
the C + solution while the solid lines represent the C_ solution. 
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We concentrate first on the travelling wave. Since the wave 
has a propagating character we interpret this mode, according 
to our model, as a stable solution, the unstable mode simply 
grows with time and it does not propagate along the thread pro- 
ducing oscillations (in the absence of longitudinal flows). The 
thread radius and the wavelength of the perturbation are known 
from observations, therefore we calculate the critical Alfven 
speed to have stable solutions. From Eq. ([33} it is found that 
v c . w 17.7 km s . This is the lower bound to have stable modes 

Ap 

whose wavelength is equal to or smaller than 3000 km. The 
Alfven speed in the thread might be higher than 17.7 km s~', 
and the exact value is calculated using the period of oscillation 
known from observations, and the dispersion relation. From the 
two possible modes, only the solution associated to the lower 
interface can be matched with the parameters determined from 
observations. Figure [6] shows the dependence of the square of 
the frequency with the Alfven speed in the thread. The Alfven 
speed that produces a match with the period of 5.4 min (and the 
corresponding phase speed of 9.5 km s _1 ) is 18.8kms _1 , and has 
been calculated using Eq. OTT i. This value corresponds to the in- 
tersection of the long dashed line, representing the reported pe- 
riod, and the thick continuous line associated to the dispersion 
relation. Note that the difference between the critical and the es- 
timated Alfven speed in the thread is rather small. 
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Fig. 8. Square of the frequency as a function of the Alfven speed 
in the thread. In this plot R = 3.3 x 10 2 km, A = 1.2 x 10 5 km 
and V° A ~ 389.1 km s _1 . The same notation as in the previous fig- 
ures has been used. The long dashed line corres ponds to the frequenc y 
(P = 4 min) derived from the observations of Okamoto et al. (2007). 
The dashed lines correspond to the C+ solution while the solid lines 
represent the C_ solution. 



200 [ 



150 - 



100 - 



50 - 



-.10" 6y-1 0" 
A (km) 



8x10' 



1x10" 



Fig. 7. Growth rate (r = 1/M) as a function of the wavelength of the 
perturbation. In this plot R = 10 2 km and v Ap = 18.8 km s _I . The dashed 
line denotes the wavelength, A c = 3.2x 10 3 km, where the mode changes 
from stable to unstable. 



Once we have an estimation of the Alfven speed in the thread 
we calculate the critical wavelength to have a stable situation. 
Using Eq. ( l32l we obtain a value of 3200 km. This value is 
slightly larger than the wavelength of the propagating wave, 
3000 km. In Figure [7] the growth rate is plotted as a function 
of the wavelength, being the lower value the critical wavelength. 
Close to the critical wavelength the growth times can be very 
large but for A > 10 4 km the growth rate is almost constant with 
a value of 20 s. Using this growth rate, the value of zo, and the 
value of the amplitude vo assumed to be of the order of 1 km s , 
we estimate from Eq. (l39l that the mixing time is t m - 48 s. 
This value is one order of magnitude smaller than the estimated 
lifetime of the thread, which is around 540 s (9 min). In this par- 
ticular case the instability is too efficient and according to our 
model this thread should disappear in less than one minute. 

Another example of thread oscillations has been investigated 
bv lOkamotoet al. (2007). Using Hinode observations these au- 
thors found transverse mot i ons of threads in an active region 
prominence. iTerradas et al.l (120081) performed a seismological 



analysis of the same event. The radius of the threads are in the 
range 180-330 km and their length is 1700- 16, 000 km. The es- 
timated minimum length of the magnetic tube is L — 10 5 km. We 
concentrate on the thread of length 16,000km since this struc- 
ture has the smallest longitudinal motion (the horizontal flow 
velocit y is around lSk ms 1 ). Contrary to the situation stud- 
ied by iLin et al.l d2007l) the oscillations seem to have a stand- 
ing nature, and the thread is moving verti cally with a period 
of 4 min. Another difference with respect to lLin et al. (2007) is 
that we have information about the maximum wavelength that 
satisfies the line-tying conditions. It corresponds to the funda- 
mental mode (A = 2L) meaning that for a fully filled thread 
A - 2 x 10 5 km. However, we c an take in to account the fact 
that the tube is only partially filled. Soler et al. (2010) calculated 
the frequency of oscillation for a thread of length L p (the total 
length of the tube is L). From the comparison of this frequency 
with the frequency of oscillation of the fully filled thread they 
found that the effective wavelength that should be used is A — 

n J(L - Lp] L p . This expression is valid as long as L v /L < 0.4. 
For the thread we are considering we have L p = 16,000 km 
while L = 10 5 km, (L p /L - 0.16), and the corresponding ef- 
fective wavelength is A - 1.2 x 10 5 km. This is a simple way 
to mimic the effect of a partially filled tube. Using this effec- 
tive value for the wavelength we calculate from Eq. (|33T > that the 
critical Alfven speed to have a stable situation is 389.1 kms~'. 
Therefore, if the Alfven speed is larger than this value the thread 
is going to be always stable with respect to the Rayleigh-Taylor 
instability. Figure [8] shows the dependence of the square of the 
frequency with the Alfven speed in the thread for the previous 
example. Again, as for the situation studied in ILin et al.l (|2007) 
only the mode associated to the lower interface can be fitted with 
the reported period, and the inferred Alfven speed in the thread 
is 525.8 km s . This speed is in the range of values estimated by 
ITerradas et al.l(l2008l) . However, the interpretation in terms of the 
interface mode does not agree with observations since for this 
particular example the whole thread shows a clear kink-like ver- 
tical motion. Again the effect of gravity seems to be too strong 
in the slab model. 
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5. Discussion and conclusions 

In this work threads have been modelled as Cartesian slabs. 
The presence of gravity acting perpendicularly to the axis of 
the structure does not allow us to have a simple equilibrium us- 
ing a more appropriate geometry. The most elemental models 
in cylindrical geo metry involve the pres e nce of magnetic twist 
(see for ex ample Lerc he&Lowi 119801: iLow&Zhangl l2004t 
iPetrie et all 120071: Blokland & Keppeni. 12011 ah . To simplify 
things and to have an analytical treatment of the problem, in- 
stead of a thread with cylindrical cross-section we have adopted 
the slab configuration. We have supposed that the slab is infinite 
in the ignorable direction, and more important, we have assumed 
that there is dense material along the whole tube since the gravity 
force prevents to have a simple analytical equilibrium between 
the dense and the evacuated part of the thread along the axial 
magnetic field. These features of our model certainly affect the 
frequencies of oscillation of the stable modes and the growth 
rates of the unstable modes. In addition, the effect of gravity is 
overestimated in the slab configuration because the gravity force 
is perpendicular to all the points on the interfaces separating the 
two fluids (see Figure [TJ. In a more realistic three-dimensional 
model, such as a cylindrical tube, the overall contribution of the 
gravity force perpendicular to the interface is weaker. This can 
be understood from Figure [9] gravity is perpendicular to the in- 
terface only at points A and B, while at point C the projection 
across the interface is zero. However, we have to bear in mind 
that a realistic geometry could be anything between a slab and a 
flattened cylinder, which might be closer to a slab of finite width 
rather than to a pure circular cross-section cylinder. 

We have shown that the mode that fits with the frequency 
reported in the observations is always associated, according to 
our model, to the lower interface and does not have a global 
nature. This is in c ontradiction with the observed motions of 
lOkamoto et all d2007l) displaying a global vertical displacement 
of the whole thread. If the effect of gravity on the eigenfre- 
quencies is smaller, as we expect for a more realistic configu- 
ration, then a match in frequency for the global mode is possible 
(the upper curve in Figure [8] might intersect the dashed line). 
Additionally, the growth rate of the unstable mode will increase, 
which is one of the questionable results of the slab model since 
the growth of the instability is too fast, according to the compar- 
ison with the observations of Lin et al.l d2007l) . 




Fig. 9. Sketch of the three-dimensional cylindrical model representing 
a prominence thread. The projection of the gravity force varies along 
the circular interface. 



We have also adopted in our analysis the incompressible as- 
sumption for the perturbations. It is known that, in general, the 
effect of compressibility helps to stabilise the unstable modes. 
Since gas pressure must change with position to have hydro- 
static equilibrium, the analytical investigation of the compress- 
ible case is much more complicated than the incompressible sit- 
uation due to the presence on non-constant coefficients in the 
equations (for example, the sound speed would change across 
the slab). However, even under the compressible assumption, in 
the regime we have considered in this work, k x <K k y , motions 
are essentially incompressible. A clear example o f this situation 
is the kink mode in the compressible slab (see Arregui et all 
2007) or t he kink mode in cylin drical geometry in the thin tube 
limit (see iGoossens et all l2009h . Hence, incompressibility is a 
reasonable approximation that has the additional advantage that 
the problem can be studied fully analytically, this has allowed 
us to derive the algebraic dispersion relation given by Eq. (l22l . 
Notice also that the choice of a step function for the density is 
convenient to avoid a situation where the Alfven speed smoothly 
changes with height, since this would imply the existence of res- 
onances for the case with ky + 0. 

From the analysis of the dispersion relation and the eigen- 
functions we have shown that, contrary to the purely magnetic 
case, the global transverse kink mode does not necessarily ex- 
ist for the whole range of parameters in the stratified model, 
even in the stable regime. This might have some important con- 
sequences, since other types of motions associated to the indi- 
vidual interfaces, or of mixed nature between surface and global 
kink waves, are feasible according to our model. Thus, from the 
observational point of view it might be interesting to try to iden- 
tify such kind of motions since their detection can give more 
clues about the real physical conditions in threads, specially 
about the interplay between magnetic and gravitational forces. 

Using the slab model we have demonstrated that the system 
allows a mode that is always stable, independently of the value 
of the Alfven speed in the thread, while there is another mode 
that is Rayleigh-Taylor unstable but becoming stable when the 
Alfven speed is increased. Hence, gravity might have a strong ef- 
fect on the modes of oscillation of horizontal threads, depending 
on the strength of the magnetic field. In t he case of t hread s in qui- 
escent filaments like the one studied bv lLin etalj (120071) . where 
the Alfven speed presumably low, too short lifetimes are found 
using the slab model, i.e., the instability is too efficie nt. On the 
contrary, in active region prominences (like in Oka moto et all 
2007), the stabilising effect of the magnetic tension might be 
enough to suppress the Rayleigh-Taylor instability for a wide 
range of wavelengths. We have also shown that the instability 
can be used to establish the minimum Alfven speed in the thread, 
an approach that is different from other seismological analy- 
sis of thread oscillations (se e for e xample Terra das et all [2008; 
iLinetall l2009t ISoler et all l2010h . Nevertheless, more elabo- 
rated models need to be explored in detail before a reliable com- 
parison of periods and thread lifetimes with observations can be 
performed. 

Finally, we have to remark that the classification of the 
modes and their spatial structure strongly depends on the slab 
configuration. As we have already mentioned, the sausage mode 
found in our system when the magnetic field is dominant does 
not exist as a trapped mode in the cylindrical tube. This mode 
is precisely becoming unstable when gravity dominates over the 
magnetic field. This means that the geometry can have an im- 
portant effect on the spectrum of the mod es, and a careful anal- 
ysis using more realistic geometries (see Blokland & Keppens, 
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for flattened cylindrical models) is required to under- 



stand well the effect of gravity in horizontal threads. 
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